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On : The solid nematic equivalent of the Fredericks transition is found to depend on a critical field 

rather than a critical voltage as in the classical case. This arises because director anchoring is 
principally to the solid rubbery matrix of the nematic gel rather than to the sample surfaces. 
O ■ Moreover, above the threshold field, we find a competition between quartic (soft) and conventional 

harmonic elasticity which dictates the director response. By including a small degree of initial 
director misorientation, the calculated field variation of optical anisotropy agrees well with the 
conoscopy measurements of Chang et al (Phys. Rev. E56, 595 (1997)) of the electro-optical response 
of nematic gels. 
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I. INTRODUCTION 

| In nematic elastomers and gels, liquid crystal ordering and the nematic director are coupled to the mechanical 
<**] i degrees of freedom. Both symmetric shear strains and antisymmetric deformations couple to the director. Since the 
macroscopic shape of the polymer network naturally depends on the chain anisotropy direction, the coupling gives 
rise to some remarkable features of the elasticity: shape can change spontaneously by well over 50% on changing 
temperature through the clearing point, shears can induce director rotation, and some shear deformations cost no 
change in energy in the ideal case. We call this last effect "soft elasticity" . It is unique to nematic networks and other 
elastic solids where a non-elastic internal degree of freedom is coupled to the mechanical degrees of freedom. The 
additional symmetries leading to softness were predicted phenomenologically by Golubovic and Lubensky and from 
statistical mechanics of large deformations by Bladon et al || . Depending on the thermo- mechanical history, nematic 
elastomers have been shown to demonstrate either extreme softness or substantial deviations from soft response ||. 
Systems with such residual resistance to deformation have been called semi-soft Q. They are qualitatively like soft 
00 ' elastomers in that the same particular modes of non-trivial deformation cost a very small elastic energy, compared to 
, other non-specihc deformations. 

Conventional nematics couple strongly to electric fields (E), the director aligning parallel to the field in the case of 
positive dielectric anisotropy Ae. When the director is anchored at boundaries and the sample is uniformly aligned 
perpendicular to the direction to be taken by the electric field, there is a transition (Fredericks) at a finite critical 
voltage, Vf, above which the director is deflected by the field. The electric Fredericks effect and its analogues are 
the basis of most liquid crystalline displays. Director rotation to lower the dielectric energy in the bulk becomes 
non-uniform if the boundary anchoring of the director is to be respected. In this case the Frank elastic energy density 
q , ii^(Vn) 2 opposes the field effect, where if is a Frank constant and n is the nematic director. For a cell of thickness 
d with strong anchoring on the boundaries the Frank energy is of order ^if (Vn) 2 d ~ per unit area of cell (where 
Vn ~ i, crudely). The corresponding electrical energy per unit area is, roughly, ^e Q AeE 2 d. These are comparable 
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for a field of Ep ~ \\l - \ e or an applied voltage Vf = Epd = \J K/e Q Ae, where the Fredericks transition occurs. 
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The response of the nematic director to applied electric fields in nematic elastomers and gels raises questions about 
the role of the elastic properties of the polymer network in limiting the director response. Likewise, the reorientation 
of the director by an external field can lead to elastic deformations. 

Nematic elastomers have been seen to respond to modest electric fields Q by changing their shape when they are 
geometrically unconstrained. At first sight this is unexpected since the energy scale for shape change in a solid is 
of order /j,, the shear modulus. Balancing this with the electrical energy, the characteristic fields required to change 
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shape should be E ~ y e ^ Ae j ~ 10 7 V/m for a typical rubber modulus /i ~ 10 5 J/m 3 and a substantial dielectric 

anisotropy Ae ~ 6. This is a large field, whereas shape changes were observed for modest electric fields. However, 
because nematic elastomers can suffer certain shape changes with little or no energy cost, which can potentially 
be achieved most easily in a completely unrestricted sample, such observations add weight to the concept of soft 
mechanical response. 

More recent experiments of Chang et al || show an electro-optical response in a nematic gel in the usual constrained 
Fredericks geometry. They observe the director response by the method of optical conoscopy, and find an elevated 
threshold field rather than voltage for director reorientation. They analyze their data in the classical Fredericks 
manner. However, to explain the threshold field, they invoke both a bulk anchoring mechanism and an intrinsic, 
material-defined length scale for gradients of the director, rather than the sample thickness, which sets the length 
scale for a simple nematic. Given that finite fields were required, Chang et al were not observing pure soft elasticity of 
the gel. It is not clear that the nematic gel can react softly in the Fredericks geometry, since the delicate combination 
of shears required to accommodate the director at no cost of elastic energy is frustrated by the cell electrodes, to 
which the nematic solid must conform. 

At this point it is important to recall another physical system, prepared in a way similar to the swollen nematic 
gels of Chang et al ||. Polymer-stabilized liquid crystals (PSLC) M are formed from a mixture of common nematic 
and polymerizable monomers. On polymerization, complete phase separation occurs and, when the polymerization 
occurs in an aligned geometry, the resulting polymer fibrillcs preserve the direction of anisotropy by creating a large 
amount of internal boundary. Thus, the phase-separated polymer mesh provides both a bulk anchoring mechanism 
and an intrinsic small (micrometer size) length scale for director response. However, the response time for director 
reorientation for both field-on and off cases is fast (which is one of the reasons for PSLC use). The compounds used 
by Chang et al were highly miscible and had a slow response dynamics, characteristic of nematic gels. We assume, 
therefore, that their experiment dealt with an aligned homogeneous gel, composed of cross-linked single molecular 
strands of nematic side chain acrylate polymer, swollen by a solvent of similar mesogenic molecules, rather than a 
PSLC system. 

In this paper we reconsider the experimental observations of Chang et al. We seek to explain the new electro-optical 
transition in terms of a homogeneous response of the director, rather than invoking a new small length scale. The 
elevated field threshold and the director response above threshold are accounted for by the coupling of the director to 
the elastic polymer network. This analysis is partially supported by some new experimental observations indicating 
that the samples respond to the electric field on a scale comparable to the sample thickness, rather than on a much 
smaller scale. To summarize this new analysis, we find that the free energy density for director rotation through angle 
9 in a nematic elastomer is schematically represented by 

1 * \ K (f ) 2 " r° AsE2 sin2 9 + ( I sin2 9 + \ sin4 e ) + ■ ■ ■ ■ w 

The Frank and the dielectric terms are as in classical nematics. The rubber-elastic /i terms are the resistance the 
gel network presents to director rotation, only the quartic term being present in the ideal soft case. The parameter 
A measures semi-softness, that is the residual harmonic elastic resistance. In elastomers and gels, the director has a 
bulk as well as surface anchoring - there is an elastic resistance even when directors are rotated uniformly (equivalent 
to a massive, finite-energy q — > normal mode). One sees that the uniform system would rotate in response to an 

(a X 1 / 2 

electric field only when the threshold E* ~ I - ^ e J is exceeded, a characteristic field rather than voltage because 

of the bulk anchoring. Additional effects arise from gradients of n in order to satisfy the boundary conditions and 
from the necessity for the sample to preserve its overall mechanical shape, but the main effect is captured by the 
above argument. We shall describe both effects in this paper, first addressing the response of a uniform system and 
then examining the role of constraints and non-uniform deformations. 

II. THE ELECTRO-EL ASTIC-NEMATIC ENERGY 

We start with a model of a macroscopically uniform sample of nematic elastomer in which the director is aligned 
parallel to the plane of a dielectric cell. When subjected to an electric field perpendicular to the cell plane (assuming 
positive anisotropy Ae), the director is forced to rotate upwards towards homeotropic alignment as in the "traditional" 
Fredericks effect. We assume for simplicity that the director rotates in the x-z plane and is parameterized by an angle 
9, which is a function of coordinate z only, 9 — 9(z), but not of x, just as in a classical Fredericks transition (see 
Figure [l]). Mechanical constraints of (i) rubber incompressibility and (ii) elastic compatibility then impose strong 



2 



limitations on the number and types of possible elastic deformations. Since the director orientation is coupled to the 
elastic network, director rotation in the x-z plane causes deformations X xx and X zz (and X yy from the incompressibility 
constraint DetX =1), and especially the shears X xz and X zx . All these components of strain are possibly functions 
of z. There are restrictions due to mechanical (in fact, geometric) compatibility [pi, which can be expressed in the 
following simple way. Since the deformation (Cauchy strain) tensor is a derivative Ay = dRi/dR° with R° and R 
being positions of a material point initially and after deformation, then one can obtain the second derivative in two 
possible ways: 



d 2 R l 
dR°dR° 



dX ik 
dR° 



Making some elements Ay functions of z induces others to be functions of x. For instance, consider the component 
of shear X zx {z): we have dX zx /dz — dX zz /dx, i.e. the extension X zz must be a function of x, which we assumed 
prohibited. Hence X zx = 0. Similarly, assuming the local extension X xx (z) leads to dX xx /dz — dX xz /dx, i.e. X zx (x), 
another prohibited rr-dcpendence. Hence X xx — const. However, incompressibility demands that X zz X xx X yy = 1. 
Therefore, the constant extensional strains in the cell plane lead to the conclusion that X zz — const as well. Analyzing 
the strain tensor, we conic to the conclusion that in a constrained cell geometry depicted in Figure [t] the only possible 
elastic deformation is the shear X xz , a possible function of z. Thence the overall strain tensor is 




Nematic rubber elasticity depends on the anisotropic step-length tensor £jj = £±6ij + (£i\ — l\_)riirij of the network's 
polymer chains. Before deformation the initial director is aligned along the cell, n x = 1, and the corresponding step- 
length tensor is called l^. When the director rotates, it can be parameterized by the angle 9(z) via n x — cos8; n z = 
sin# and the corresponding step-length tensor £ is a function of 9. It is useful to introduce a ratio r = £\\/£j_, the 
measure of chain anisotropy. Main-chain liquid crystal polymers may have r > 10. The siloxane side-chain polymers 
used by H. Finkclmann |3] have r ~ 2 — 3. The acrylates used in the experiments of G.R. Mitchell Q and of Chang 
et al || have a much lower anisotropy of the backbone, r < 1.3 — 1.5. 

The cell boundaries impose additional nematic and mechanical constraints: (i) as in the classical Fredericks effect, 
the director should respect the anchoring conditions, 8(0) = 9(d) — and (ii) no overall macroscopic displacement 
along x is possible, so the strain X xz (z) must have an integer number of oscillations (see the right of Figure |l|; this 
favors a full wavelength Fredericks transition, in contrast with the half- wavelength Fredericks transition of the ordinary 
nematic, left in the picture). 





xz 




X 



FIG. 1. Field-induced director rotation in a conventional nematic (left) and nematic elastomer (middle) anchored at the 
surfaces z = and z — d, with an electric field E applied across the cell. The shear strain X xz accompanying the director 
rotation in nematic elastomers is shown on the right. The conventional Fredericks effect has one half wavelength of director 
rotation between the plates, while the solid nematic Fredericks effect has the full wavelength. 

With all these preliminary restrictions and conditions, the full free energy density of the system takes the form: 



/ - \, Tr 



-H A (sin 6 - X xz cos 9) 2 - ^e AeE 2 sin 2 + ^K(Vn) 2 
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The first term is the ideal nematic rubber elastic energy density [Q and can lead to the soft elastic response. The 
second term is the appropriate non-ideal, "semi-soft" contribution observed in many elastomers and, in particular, 
in those formed in a field-aligned monodomain state p0[ . In elastomers /i gives the rubber energy scale, essentially 
a shear modulus at small deformations. The dimensionless parameter A measures the (usually small) semi-softness. 
Inserting the deformation A and the director rotation 9 (through £) , these two rubber- nematic terms take the form 
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-\ 2 XZ + A(sm9 - X xz cos6») 2 



(3) 



where we henceforth absorb the constant 3/i/2, the energy of the relaxed state, into /. The last two terms in equation 
(^f) are conventional for liquid crystals. We adopt a simplified form of the dielectric term, with the field E rather than 
the displacement D. The reasons for and the applicability of this simplification are discussed in the Appendix. 



III. ANALYSIS OF UNIFORM DIRECTOR ROTATION 



Without any assumptions of small deformations or angles, the optimal shear strain is given by the minimization of 
the free energy density f ru b(9,\), Eq. (|^), with respect to X xz at a given 9. One obtains, after some straightforward 
algebra: 

, , _ (r- A- 1) sin 2(9 

2[l + Ar+(r- A- l)sm 2 9Y W 

where the local director rotation angle is, in fact, a function of z, 6 = 9(z). Substituting this back into the free energy 
density gives an effective free energy density, /*, depending only on the 'liquid crystal' variable 9. 

rW = ^(-(Vn) 2 -f^sin^+ (5) 
(r- l) 2 sin 4 (9 A r 2 sin 2 9 



1 + (r - 1) sin 2 9 (1 + (r - 1) sin 2 0)(1 + Ar + (r - A- 1) sin 2 6 

The effect of the underlying rubbery network is expressed by the term in square brackets, f* ub (9), an effective penalty 
for the uniform director rotation imposed by the elastic gel network. This expression is arranged in such a way that 
the contribution from the ideal soft nematic rubber elasticity is separated from the additional semi-soft part, the term 
proportional to the parameter A. 

One can notice that at small deformations, 9^1, only the semi-soft correction contributes to the elastic energy, 
the first term in brackets being proportional to the higher power 9 A (appropriately expressing the softness of the 
ideal rubber-nematic response). Therefore, the main threshold for the uniform director rotation is determined by the 
semi-soft parameter A, which controls the counter-torque to the electric field. Note, that any effect of the field is 

1/2 

scaled by the rubber modulus factor, yielding a reduced dimensionless electric field £ = (eoAeE 2 / /i) . In many 
cases this makes the response extremely small. However, the highly swollen network used in the experiments (up 
to 90% solvent) may well have a much lower modulus (as will be indeed confirmed in Section |vi|) . Also, as we shall 
see below, the scaling depends on the important polymer anisotropy (r — 1) too, so that the materials with small 
backbone chain anisotropy have, in fact, a stronger response to electric fields. 

In spite of the obvious non-linearity of the effective free energy density (||), it is possible to solve exactly for the 
optimal director angle 9(E) by minimizing f*(9). The uniform, field- induced rotation of the director is 



•2/1 Ar + 1 , ^Ar + 1 

sm 9 = ■ ; 6 

r-Ar-1 y/r-1- £ 2 y/r ~Ar-l v ' 

or ^ 2 ^ (r _f' gc2)2 [g/gc-l], (7) 

where the second expression shows the variation of angle just above the bulk threshold. The appropriate reduced 
variable x = £/£ c — 1 is used near the threshold field, £ c . The value of this threshold is: 

Ar 2 



4 



and is determined by both the semi-softness parameter (i.e. the director locking in the bulk) and the polymer chain 
anisotropy r. The solution Eq.(||) reaches its upper-bound at the field value 

fig= ( r -l)»( r + l)_ A r-2-A- > (Q) 

when sin (9 = 1 and the director rotation is complete. All this is in contrast with a conventional Fredericks effect, 
where there is no barrier for the uniform bulk director rotation and, if no boundary anchoring were involved, the 
director would switch to 9 = ir/2 immediately at E ^ 0. 

The conoscopy technique, used in experiments [||, reveals the change in optical path length for extraordinary rays 
traversing the sample 8 = j(n e — n e ff)i with d the cell thickness, and A the wavelength of light. The effective 
refractive index depends on the current director orientation 9(E) (we assume uniform rotation in this Section) 

<ff = 2 + , f' n K ■ (10) 

n o + ( n e ~ n o) Sm " 

2 _ 2 

Near the threshold, when 9 <C 1, the change in optical path is simply 8(E) ~ 4 n e 9 2 (£). A plot of 8(£) for 
d/X ~ 200, n e = 1.74, n Q = 1.54, very small semisoftness A = 0.001 and increasing values of the chain anisotropy 
r = 1.1, 1.5 and 2.5 is given in Figure 2(a). Chains with higher anisotropy r show a stronger resistance to an external 
torque, as one can see from equation (||). At fixed parameter A, the director response to the field is stronger when the 
chain anisotropy r is reduced, as can be seen from the (expanded) equation (H). However the underlying resistance 
to rotation is because of the semi-soft, quadratic (A) terms in the rubber-nematic free energy, which is why the bulk 
threshold scales like £ 2 ~ A. Only at larger rotations do the quartic terms (of order 1 rather than A) dominate, 
accounting for the scaling of the saturation field E\. Figure 2(b) shows the optical path 8(£) at fixed chain anisotropy 
r = 1.5 and increasing degree of network semi-softness, A = 0.001, 0.02, 0.05, 0.1. At small angles 8 cx 9 2 and 
hence 8 ~ (£ — £ c ) [see Eq.d?])]. Therefore increasing A (and thus £ c ) increases the linear slope of the plots after the 
threshold, as seen in Figure 2 (b). 




FIG. 2. Change of extraordinary ray path length 8(8) versus electric field, (a) Small semi-softness A = 0.001 with 
increasing chain anisotropy r — 1.1, 1.5, 2.5 giving less rapid growth of S with £. (b) At fixed anisotropy r — 1.5 with increasing 
A = 0.001, 0.02, 0.05, 0.1 showing an increase in the threshold field £ c with increasing A. 

Figure 3 uses the further reduced field x = £/£ c — 1, and the optical pathlength scaled by (r — 1 — £ 2 ) 2 /£ 2 , as 
suggested by equation (fjj). Deviations from the initial linear response are evident for high threshold fields. From 
equation (||) for sin 2 9 we see there is a qualitative change in the field response for £ 2 = r — 1, that is where the 
semi-softness A = (r — l)/r. This condition is met either at low fields in networks of low anisotropy, r — 1 small, or 
for fairly strong semi-softness, A large. At £ 2 = r — 1 we have £ 2 — £ 2 . There is no response to a field until it reaches 
£ = £ c , whereupon there is immediate, discontinuous switching 9 — * 7r/2. 

Clearly, the existence of a threshold in field, rather than voltage, and quite a steep increase of the director rotation 
[and consequently of the observed optical path 8(£)] are described well. However, there are two important aspects still 
missing: (1) the analysis of boundary conditions and non-uniform deformations, which leads to a Fredericks threshold 
voltage in conventional nematics, and (2) the experimentally observed small increase in rotation before the principal 
threshold, a "precursor foot" which we shall associate with a small degree of random misorientation in the initial 
director distribution. 
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FIG. 3. The conoscopic response for the conditions of figure 2(b). The field is further reduced to x — £/£ c — 1, and optical 
pathlength is scaled by (r — 1 — £% ) 2 /£c , which collapses the plots of Fig. 2b for various A onto the single curve, Eq.(Q), in the 
region of the transition, i<1. 



IV. CRITICAL ANALYSIS OF THE TRANSITION REGION 



The analysis of the previous Section assumes that there is no relevant spatial variation in either 6{z), or X xz (z): 
the solution for 9(E), equations (^) and (|t]), has been obtained without taking into account the director-gradient 
Frank elasticity. This is a reasonable assumption at fields high above the threshold, where the texture is coarsened 
and the bulk is capable of achieving its optimal value of 9, equation (^|), which we shall call 9q (similarly, for the 
conventional Fredericks effect do — > tt/2 and the profile 9(z) coarsens at high fields). Near the threshold one has to 
be more careful and examine the effect of boundary conditions (which is the only existing barrier for the conventional 
Fredericks effect). 

Near the threshold we can safely assume 8 < 1 and the free energy density takes the form, standard for elliptic- 
function analysis (see Appendix) 

lr»\m'-\pfi + y, (id 



with the coordinate z scaled by the natural length scale £ = y/K/fi. For a typical rubber modulus, as in experiments 
pfj| , £ ~ 10~ 8 m. Parameters p and q are obtained by expanding the full free energy density (||): 

A r 2 

x9 Ar 2 ( r - 1 - A 

- 1) r- \r - H ; — 

v ; 1 + Ar V 1 + Ar 

In q we have put £ 2 ~ Ar 2 / (1 + Ar); that is we are close to the threshold, and have neglected other terms discussed 
in the Appendix. The details of this type of analysis are also given in the Appendix of paper |fTo|| : the director 
cannot achieve its local optimal orientation #o = yjp/q, being held by boundary conditions. Instead, the director is 
only having a modulation of amplitude m (< 6q) across the cell. The resulting "order parameter" rj(E) = 6 m /8o 
increases from zero towards 1 as the texture coarsens. The number of periods of the elliptic function that describes 
the variation 9(z) is 4 J^^ K \Jv^ V 2 ) with K the complete elliptic integral of argument rf/(2 — rf ). The minimal 

number of periods in our case is one (in contrast with the usual symmetric Fredericks effect, where it is a half ) , so 
one obtains the director rotation angle 9 m as the solution i](E) of 

1 -m = 4€ as) 



1 d U*J«¥e* 



4^2 VK\ ^l + Ar y 2V2V K 

The expression in parentheses is the corresponding right-hand side for the conventional Fredericks effect, with no 
underlying rubber elasticity fi and A = (and with the half- rather than full-period of director oscillation across the 
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cell thickness); it depends on a voltage V = E d. The real threshold is shifted from that of the uniform analysis of 
last Section, equation (||), by the additional contribution from the domain walls (which is the only contribution in 
the conventional Fredericks effect): the solution of ( p"3| ) first appears at the minimal value of K = ir/2 at 77 = 0, i.e. 
p* = (2</d) 2 : 

_ 17 fi \ A~ /27r\ 2 K 
~ V \s As J 1 + Ar + \~d J e Ae 

One should note that the analysis of the threshold itself is elementary and does not require the investigation of the 
elliptic function. In the usual way, we could introduce a variation of wave vector A; (= 2ir/d) whereupon the expanded 
free energy density ( p"l| ) is / ~ ^(k 2 — p)0 2 + t<70 4 - Modulation starts when the # 2 -coefficient becomes negative for 
the first time, whence the threshold condition ([14]) is recovered. 

Returning to the elliptic analysis, near the threshold the director rotation is given by tj 2 — | (p — p*) jp* , describing 
the full-wavelength modulation with the amplitude 9 m = rj(E)9o(E). The effect of this elliptic analysis near the 
threshold rapidly disappears as the modulation coarsens, the optimal value 8q of the director rotation being achieved 
in most of the sample and the 6(E) variation becoming that of equation (|l0|). 

We said that one normally expects ([id 2 jK) ^> 1, the thickness-dependent term thus being irrelevant in the 
expression for the threshold. As discussed below, effects observed in thin samples suggest that the effective rubber 
modulus, \xA, is low in the studied material. Although this may in part be because fi is small due to dilution by 
solvent, it clearly indicates that A is small too. If the effect were simply due to the smallness of /i alone, then the 
quartic terms would also be weakened and the saturation field £ u would also be accessible. Because A is evidently 
small, it may be possible to see the remaining d-dependence of the apparent threshold in experiment and compare it 
with (Q. 



normally we expect 



Afid 2 
K 



> 1 



(14) 



V. OBLIQUE INITIAL DIRECTOR 

There seems to be no reason to assume a systematic obliqueness (pre-tilt) in the initial director. However, the 
preparation of networks, especially in thick cells, results in a residual polydomain texture, as reported by Chang et 
al. and is a well-established phenomenon in nematic gels. Even after orientation by strong magnetic fields one should 
expect a small degree of remaining random disorder, see jn|. Assuming the degree of this misalignment is small, it is 
straightforward to make provisions for it and obtain the resulting quenched averages of observable parameters, such 
as the director angle 9(£) and the optical path 6(£). 




0.1 0.2 0.3 0.4 0.5 
Reduced field S 



FIG. 4. Conoscopic path length 8 vs. reduced electric field for backbone chain anisotropy r = 1.5, semi-softness parameter 
A = 0.05 and increasing r.m.s. misalignment angle a — 0, 0.01, 0.03, 0.05, 0.1 radians. 

We do not go into details of this calculation, which involves some heavy algebra. Instead, we plot the results for 
several values of parameters. The plots in Figure [l| show the optical path 6(£) for selected values r — 1.5, A = 0.05 
and increasing r.m.s. misalignment a = 0, 0.01, 0.03, 0.05, 0.1 (the last value corresponds to a = 5.7 deg). 
The misalignment not only creates a precursor foot before the main threshold, but also changes the slope so that the 
extrapolated value of £ appears slightly reduced from E* in (^) . We emphasize that theory and physical expectations 
suggest the obliqueness a should strongly depend on cell thickness and on the magnetic field used during crosslinking. 
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VI. EXPERIMENTAL OBSERVATIONS 



We now apply the model of the previous sections to the experimental observations of Chang et al. , in light of 
further experiments performed recently. The analysis of data performed by Chang et al. includes a bulk anchoring 
effect similar to that proposed here, which produces a threshold field, rather than voltage. However, they also appealed 
to an intrinsic material length scale on the order of micrometers, which further increases the observed threshold field 
and also sets the scale for director gradients, which limit the director rotation above threshold. They reported 
observations of micrometer scale speckles in the sample after polymerization in support of this proposed length scale, 
but did not directly observe a periodic pattern of director rotation at the micrometer scale in response to the applied 
electric field. 

We prepared a series of samples of different thicknesses according the the methods of Chang et al., and performed 
microscopic observations during application of the field. In thick samples, comparable to those of we observed the 
same general elements of response that they reported, but in addition, we were able to observe small scale movements 
in the sample, particularly displacements of microscopic particles and speckles in the sample. We observed lateral 
displacements in the plane of the sample on the order of a fraction of the sample thickness, and saw no evidence of 
micrometer scale deformations in response to the field. Our observations are consistent with the deformation sketched 
on the right in figure |], which involves a combination of shear strains resulting in lateral displacement of material 
in the midplane of the sample. The deformation proposed on the right in Fig. [I] also satisfies the preservation of 
symmetry in the conoscopic image reported by Chang et al, in contrast to the rotation of the optical axis that would 
accompany the classical half- wavelength Fredericks transition. Thus, the field response proposed here is completely 
consistent with the general observations reported by Chang et al. |^|. Moreover, as described in the previous Section, 
the pretransitional response, or rounding of the transition, observed by Chang et al. can be similarly explained by a 
small random pre-tilt of the director. 
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FIG. 5. Theoretical curve of 5(E), the conoscopic length, against electric field for the present theory (solid line) for thick 
samples (uniform director rotation) to the data of Chang et al. |J . For the 62 /im thick sample (crosses) we have doubled the 
data values for the optical pathlength 8 to collapse the data from two samples onto the same plot. For the 125 /im thick sample 
(squares) the data are unsealed. This is the same scaling used in Fig. 1 of M. 



Taking the published data of Chang et al., and using the model of Section [II appropriate for thick samples, along 
with pretilt of the director, we present a fit of the present theory to the data in figure ||. We plot our theoretical 
prediction of S using 9 from ^| and taking the values r = 1.95, A = 0.06, n = 1.6 x 10 2 J/m 3 , and a — 1.7 degrees 
for the network, in addition to the sample characteristics specified by Chang et al.. We see that this theory can 
reproduce the experimental observations with reasonable values of sample thickness etc. fl3|| . The values of semi- 
softness parameter A and of the random local pretilt angle a are sufficiently small not to contradict the known facts 
about liquid crystalline gels. The value of backbone anisotropy r, slightly higher than that reported for solvent free 
polyacrylate elastomers (~ 1.44 p|Jl4|), is likely due to the large amount of nematic solvent present in the gel and 
increasing the nematic field. The solvent (up to 90%) can also explain the very small value of rubber modulus fi. In 
fact, it is mostly due to this weakness of swollen elastic network that we can observe a noticeable director rotation, the 

1/2 

angle of which depends on a ratio E (e Ae//i) . However, in spite of this plausibility, independent measurements of 
these parameters are needed to test the presented theory more explicitly. 

The value of the gel modulus \x, used in the fitting above, is substantially lower than that expected for ordinary 
rubbers. It necessitates re-examination of the characteristic length scales in our system. The natural nematic elastomer 
scale £ = y/ K//J, (see Eq. (pi]) and, more comprehensively, the review |Q) is now of the order 10~ 6 m. This scale gives 
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the characteristic thickness of domain walls, two of which are present in the full-wavelength Fredericks geometry, 
Fig. 1. Hence one expects that thin samples, with d only of order of several micrometers, should deviate from the 
description proposed here: In thick samples, our electro-mechanical response is driven by a reduction in electrical 
energy because of director rotation toward the field with very little elastic penalty, since the deformations are soft 
in the uniform regions above and below the midplane. These regions are separated by a domain wall in the middle, 
where the shear strain reverses sign. When the sample thickness is reduced, this wall and the two non-uniform regions 
near the cell boundaries would start merging. This raises the elastic energy of the full-wavelength (nematic elastomer- 
like) Frederiks transition above that of the half-wavelength (liquid nematic-like) one. Indeed, in thin samples (8 to 
25 micrometers in thickness), we observed the different response, almost identical to the ordinary half- wavelength 
Fredericks transition. This was evident in the conoscopic observations, by the rotation of the optical axis in response 
to the applied field, as in low-molar mass nematic liquid crystals. 

To satisfy the boundary conditions with this fully non-uniform director rotation, the predominant shear strain in 
the bulk of the sample must be zero. The energy penalty and the counter-torque for the action of electric field are then 
provided by Eq.(||) with X xz = 0. This is the regime of "couples without strains", with a bulk barrier for the director 
rotation, first proposed by de Gennes Jl^l . More detailed experimental measurements and theoretical analysis will be 
carried out to explore this kind of transition in thin samples, and to determine the thickness at which the crossover 
from thin to thick behavior occurs. This will provide a detailed understanding of the combined director and shear 
strain structure within the boundary layers and give an independent measure of the effective rubber modulus /i via 
the characteristic length scale £. 

Another aspect that needs to be re-examined for a system with very low modulus [i is the threshold condition, 
Eq. (|l4|) . In order to predict the field threshold E*, rather than the voltage V* — E.d, we need to maintain the 
inequality Afid 2 /K 1. For the values used in fitting the experimental data in Fig. 5 (for the smaller d ~ 60 /im) 
we still obtain a ratio of order 3 x 10 3 and the conclusions about the threshold remain valid. However, this ratio will 
become of order unity for d ~ 10~ 6 m, where the ordinary voltage threshold should again become relevant. 

VII. CONCLUSION 

In summary, we propose a model for the Fredericks transition in nematic polymer gels. Especially for thick samples, 
the polymer network serves both to set the threshold field for the transition, and to control the director response 
above threshold. It effectively replaces the role of the sample surfaces in the ordinary nematic Fredericks transition. 
We can reasonably explain the experimental results of Chang et al. without appealing to a new internal length scale 
for director gradients, and this element of our explanation is supported by qualitative experimental observations of 
macroscopic lateral displacements in response to the applied field. 

Because two coupled physical fields, the director angle 9 and the shear strain X xz , are at play, it is possible that 
there are two different characteristic timescales for relaxation. Depending upon the disparity of these, one could at 
short timescales possibly see another decay regime before the final, slowest exponential decay reported by Chang et 
al.. This obtains naturally from the numerical solution of the dynamical equations in 9 and X xz resulting from Eq.(|^) 
(see®. 

We have begun to examine the remaining role of director and shear strain gradients in the sample, and our 
experimental observations suggest that this is important for understanding the transition in thin samples. Much 
experimental work remains to be done to test various aspects of our model. Finally, it is interesting to compare 
these predictions for nematic gels with the behavior of polymer-stabilized liquid crystals. The key difference is the 
permanently fixed macroscopic polymer skeleton of the PSLC system, in contrast with the mobile molecular scale 
network of a gel. The detailed analysis of the threshold, director evolution at higher fields and the role of alignment 
conditions at preparation should distinguish between the behavior of these two systems. The dynamic behavior is 
especially different, with gels responding very slowly, compared to PSLC's. 
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APPENDIX A: ANALYSIS OF DIELECTRIC PROBLEM 

In the classical electric Fredericks problem the anisotropic dielectric fluid is contained between two parallel plates kept 
at a relative voltage V The analysis can easily be generalized to the elastic-electric case. One needs to use the 







electric displacement vector D = (0, 0, D), which is a constant between the plates given no space charge (div D = 0). 
The field E and the displacement D are related by D = e e-E. As all nematic uniaxial tensors, the matrix of dielectric 
constants is expressed as e# = e±5ij + (ey — e±)mnj, where the nematic director n is the principal axis. The voltage 
is, therefore, given by the integration J E z dz: 
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e (e_L + Aesin 2 9) 
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where e zz = e± + Aesin 2 9. The correct free energy, now per unit area of plates, does not derive from equation (j^), 
but from the dielectric energy density — (D • E): 
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where ng(s) is the rubber-nematic part of free energy density, given by equation 
giving the 6(z) minimizing this free energy is: 
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(A2) 



The Euler-Lagrange equation 



(A3) 



where the constant D has to be determined from equation (|Al| ) and is thus a functi ona l of 9{z). The boundary 
conditions are 9' — at 9 = 9 m if there is a modulation with amplitude 8 m . Then Eq. (A3) is easily integrated once 
to give: 



D 2 

5(sin 2 9 m )) + 



1 



1 



ej_ + Ae sin 2 9 ej_ + Aei 



(A4) 



where u = z/£ is the length reduced by the rubber-nematic correlation length £ = yl Kj [i. The full problem is solved 
by further integrating over 1/4 period and yields: 
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(A5) 



where [. . .] is the right hand side of Eq. ( A4). The right hand side of Eq. ( A5) is a quarter period of the modulation. 
Integrating Eq (A4) to a general 9 at a general z, rather than to 9 m and d/4 as in Eq. (|A5|), gives 9(z) as an elliptic 
function, which then finally determines D via Eq. (Al) and thereby in Eq. (A5) itself. 

At the vicinity of the transition, all this analysis significantly simplifies since sin 2 9 < sin 2 9 m < sin 2 9 a « 1. D 

D 2 1 



and V are trivially related as D = e e±V/d. The potential g{s) 



E £j_+ Aesin 2 



driving Eq. (A4) simplifies to 



— \"p8 2 + jq9 4 with p and q given in the equation (p^). The elliptic analysis then follows as if from the simplified free 
energy Eq. (§)• 

In the uniform case considered in Section |l| the distinction between D and E is easily handled too since Eq. ( Al ) 
collapses to the obvious D = e Q (e± + Aesin 9)V/d and the electrical part of equation (A2) gives the ^-dependent 
part of the electrical energy in equation (|^) (a constant part having been suppressed). 
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